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ABSTRACT 


The Telescope Array Collaboration has observed an excess of events with E > 1019-4 eV in the data 
which is centered at (RA, dec) = (19°, 35°). This is near the center of the Perseus-Pisces supercluster 
(PPSC). The PPSC is about 70 Mpc distant and is the closest supercluster in the Northern Hemisphere 
(other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with 20°- 
radius circles indicates an excess in the arrival direction of events with a local significance of about 4 
standard deviations. The probability of having such excess close to the PPSC by chance is estimated 


to be 3.5 standard deviations. 
supercluster. 


This result indicates that a cosmic ray source likely exists in that 


Keywords: Particle astrophysics (96), Ultra-high-energy cosmic radiation (1733), Cosmic rays (329), 
Cosmic ray astronomy (324), Large-scale structure of the universe (902), Cosmic ray sources 


(328) 


1. INTRODUCTION 


Ultrahigh energy cosmic rays (UHECR) are energetic 
particles originating from outer space, having energies 
greater than 1018 eV, that impinge on the Earth’s atmo- 
sphere. Their sources remain one of the most important 
questions to be answered. The main areas of UHECR 
study are their energy spectrum, mass composition, and 
searches for anisotropy, which are expected to shed light 
on the question of sources. The Telescope Array (TA) 
hotspot (Abbasi et al. 2014) which is about 20° in ra- 
dius and centered at (146.7°, 43.2°) is one indication of 
anisotropy. Other previous signs are the Auger dipole 
(Aab et al. 2017), most of the evidence for which is in 
the southern hemisphere, and the Auger excess (Aab 
et al. 2015) near the location of Centaurus A, the closest 
active galactic nucleus. The nearby Virgo cluster, dom- 
inated by M87, does not appear as an excess in UHECR 
data. 


* Presently at: University of Californa - Santa Cruz and Flatiron 
Institute, Simons Foundation 
+ Deceased 


The TA collaboration hotspot, reported in 2014 (Ab- 
basi et al. 2014), was found in the arrival direction 
of UHECR events with energy greater than 57 EeV 
in the northern sky. In that publication, an oversam- 
pling analysis using the intermediate angular scale of 
20°-radius circles was conducted for the first 5 years of 
data collected with the TA scintillator surface detector 
array. The maximum excess appeared at the position 
(146.7°, 43.2°) in equatorial coordinates with a Li-Ma 
significance of 5.10. The probability of having such an 
excess by chance was estimated via Monte Carlo (MC) 
studies to be 3.40. However, there are no known promi- 
nent sources aligned with the hotspot within the cosmic 
ray events’ horizon. There have been correlation studies 
between the hotspot and possible sources (He et al. 2016; 
Kim et al. 2019); however, the source of the hotspot re- 
mains inconclusive. The hotspot persists in TA data 
with over 30 significance (Kim et al. 2021). 

In the Northern Hemisphere, the next closest com- 
ponent of the local large-scale structure is the Perseus- 
Pisces supercluster (PPSC). It stretches across the sky 
from the Perseus cluster to the Pegasus cluster and con- 
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tains about 16 galactic clusters and groups containing 
tens of thousands of galaxies. 

In this paper, we report a new excess of events at 
slightly lower energies than the original hotspot. While 
studying the spectrum mismatch above 10!%-> eV seen in 
the TA and Auger data (Abbasi et al. 2021a), we made 
sky maps of events for energy ranges, E > 10/94 eV, 
E > 10 eV, and E > 10/9 eV. The maximum Li- 
Ma significance appeared at the location of the PPSC 
for each energy range. We now present these three data 
sets. 

The structure of this paper is as follows. In Section 
2 we describe the TA experiment and the data used in 
this study. Section 3 has three parts. First, we present 
the Li-Ma oversampling analysis using 20°-radius circles 
to determine the significance of the excess of events in 
arrival directions. Then we show sky maps made with 
the representative elements of the PPSC. Finally, we 
describe the MC simulation used to estimate the chance 
probability of having an excess close by the prominent 
local large-scale structures. In Section 4 we summarize 
our findings. 


2. TA EXPERIMENT AND THE DATA 


The TA is the largest observatory for UHECRs in 
the Northern Hemisphere. The observatory is approx- 
imately 1400 m above sea level, and it is centered at 
39.3°N and 112.9°W, in the west desert of Utah, USA. It 
consists of a surface detector (SD) array and three fluo- 
rescence detector (FD) stations viewing the sky over the 
array. It is designed for observing extensive air showers 
induced by UHECRs using a hybrid technique. 

The SD array consists of 507 plastic scintillation de- 
tectors deployed on a square grid with 1.2 km spacing, 
covering a total area of ~700 km? (Abu-Zayyad et al. 
2012). They measure the shower footprint and lateral 
distribution at the Earth’s surface. Three FD stations, 
instrumented with 38 telescopes, are situated at the 
apices of a triangle with each station having a field of 
view overlooking the area of the SD array (Tokuno et al. 
2012). The FDs are suitable for measuring the longitu- 
dinal shower development of events to estimate the mass 
of the primary particle; however, their data are limited 
to ~10% duty cycle since they operate only on moonless 
nights. In this study, we focus on the data recorded by 
the SD array to take full advantage of its high duty cycle 
(greater than 95%). 

For this work, we used the data collected between May 
11 of 2008 and May 10 of 2019—a data-taking period of 
11 years. The event selection criteria that apply to the 
data are as follows: 

1. Energy > 1019-4 eV, 


2. Zenith angle of arrival direction < 55°, 

3. At least five SDs triggered, 

4. Shower geometry and lateral distribution function 
fit x?/dof < 4, 

5. Reconstructed pointing direction error < 5°, 

6. The fractional uncertainty of the energy estimator 
S300 < 25%, 

7. The largest signal counter is surrounded by four 
working counters: there must be at least one working 
counter to the left, right, down, up on the grid of counter 
with the largest signal. These counters do not have to 
be immediate neighbors of the largest signal counter. 

These are our standard selection criteria that have 
been used previously for anisotropy studies. The energy 
of reconstructed events is determined by the SD Sgog 
which is then renormalized by 1/1.27 as previously de- 
termined to match the SD scale to the calorimetrically 
determined FD energy scale (Abu-Zayyad et al. 2013). 
A total of 864 events meet these selection criteria. The 
energy and angular resolution of events are from 10% to 
20% and 1.0° to 1.5°, respectively (Abbasi et al. 2018). 


3. LI-MA OVERSAMPLING ANALYSIS 


In conducting the spectrum study where this work 
originated (Abbasi et al. 2021a), part of our investiga- 
tion was whether the hotspot extends down into lower 
declinations at slightly lower energies. Therefore, we 
adopted 20° angular windows for oversampling to be 
consistent with the original hotspot analysis method 
(Abbasi et al. 2014). As a result, we found new excesses 
of events in the distribution of arrival directions. 

The TA observes from 90° to —15.7° in declination 
and 0° to 360° in right ascension. To estimate the back- 
ground event rate, 10° events were generated within that 
field of view, assuming an isotropic flux and taking into 
account the geometrical exposure of the SD. The sta- 
tistical significance of the excess of the data was com- 
pared to the isotropic events at each grid point, using 
0.1° x 0.1° grid spacing. The significance was calculated 
utilizing the Li-Ma method (Li & Ma 1983): 


B (140)Non (1+a)No 1/2 
Sum = V2 [Non In (RERNE ) + Now In (SER) 
(1) 


where a = Nsim,circle/ (Nsim, total = Nim aie Here, 
Nsim,total is the total number of isotropic background 
events we generated, Ngim,total = 10°, Neim,circle Means 
the number of isotropic events inside the 20°-radius cir- 
cle at each grid point, Non represents the number of 
observed events inside the circle at the same grid point, 
and Nog = Naata,total Non, where Naata,total is the total 
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number of data in each data set. The number of back- 
ground events inside the 20°-radius circle is estimated 
by the exposure ratio a and the data, Nbg = a Nor. 

When compared to the isotropic background events, 
new Li-Ma significance results are calculated for each 
energy cut and are summarized as follows. For E > 
10194 eV, the maximum Li-Ma significance appears 
at (17.4°,36.0°) with 4.40. There are 85 out of 
864 observed events within the central 20°-radius cir- 
cle oversampling, whereas 49.5 events are expected 
as background from the isotropic assumption. For 
E > 10'%° eV, the maximum Li-Ma significance ap- 
pears at (19.0°,35.1°) with 4.20. There are 59 out of 
558 observed events within the central 20°-radius cir- 
cle oversampling, whereas 31.5 events are expected as 
background from the isotropic assumption. For Ek > 
1019 eV, the maximum Li-Ma significance appears at 
(19.7°, 34.6°) with 4.00. There are 39 out of 335 ob- 
served events within the central 20°-radius circle over- 
sampling, whereas 18.6 events are expected as back- 
ground from the isotropic assumption. 

Figure 1 shows the excesses of events with E > 
1019-4 eV, 10195 eV, 10196 eV, and 57 EeV, respectively. 
The color scheme represents the Li-Ma significance ex- 
plained above. On each map, we indicated the position 
of maximum Li-Ma significance, excess with respect to 
isotropy, with a black diamond. For the three lower en- 
ergies, the new excesses of events appear consistently in 
the direction to the left of the center in the sky maps in 
Figure 1. For the highest energy set, with E > 57 EeV, 
the maximum significance of the sky map shifts to the 
region in the figure which is the hotspot reported in 2014 
(Abbasi et al. 2014). A smaller excess remains visible in 
the direction of the maximum at lower energies. 

Next, we investigated the time variation of the excess 
by dividing the data into two time periods. As an ex- 
ample, the sky map is shown in Figure 2 for the data 
set with Æ > 109° eV. The map is shown for the first 
5-years of data adjacent to the map for the last 6-years 
of data. There is no apparent difference in the excess 
between the maps of the two data sets. Both maps have 
similar local significances, around 30, toward the new 
excess region. The distributions of the other two en- 
ergy cuts demonstrate similar local significances. This 
is indicative of a steady state excess in this region. 

The new excess of events appears in the region of the 
Perseus-Pisces supercluster (PPSC), which is one of the 
notable structures within TA’s field of view. This is 
a possible source for the observed excess. The PPSC 
is the closest supercluster other than the local super- 
cluster which we reside in. It is known that the PPSC 
has a gigantic filamentary structure, stretching for over 


300 Mpc/h (Batuski & Burns 1985). According to 
Haynes and Giovanelli (Haynes & Giovanelli 1986), the 
major portion of the foreground between the Earth and 
the PPSC as well as the space beyond the PPSC in the 
same direction, are nearly empty. Courtois et al. pro- 
vide density contour maps after making corrections for 
catalog incompleteness (Courtois et al. 2013). Figure 8 
in that paper shows the PPSC appears as an elongated 
structure of galaxies. No prominent structures between 
the Earth and the PPSC, or beyond the PPSC, are seen 
in the direction of the new excess of events in that fig- 
ure. If the new excess of events in arrival direction dis- 
tribution has an astrophysical origin, the PPSC could 
be responsible for it. 

Figure 3 shows the new excess in a series of expanded 
sky maps with the PPSC overlaid upon it. The data 
has been plotted for the three energy cuts using the Li- 
Ma significances. To show the elements of the PPSC 
on the sky, we adopt the list that includes the major 
clusters of galaxies and the major groups of galaxies that 
comprise the PPSC in reference (PPSC 2021). These are 
marked with asterisks in the three maps. The excess is 
coincident with the overall distribution of the PPSC. 

It is suggestive that the excess in the data falls on 
top of the PPSC. To quantify how often this happens 
by chance, we generate many isotropic MC event sets 
thrown according to the acceptance of the TA SD. Each 
MC set contains the same number of events as the data. 
We count as successes the number of MC sets where the 
maximum Li-Ma significance is at least as significant in 
the MC set as in the data and which occurs at least as 
close to the PPSC as the data. 

We begin by defining the directional center of the 
PPSC using the mean values of the right ascensions 
and declinations of its representative elements. It is 
calculated to be (20.9°, 27.9°) in equatorial coordinates, 
which is indicated by the blue square in Figure 3. In ad- 
dition, the center position of the excess, where the max- 
imum Li-Ma significance is estimated, is represented by 
the cyan diamond in the figure. Next, we calculate the 
angular distances between the center of the PPSC and 
the positions of the maximum excesses— (17.4°,36.0°), 
(19.0°, 35.1°), and (19.7°, 34.6°) for each energy thresh- 
old. Their angular separation is 8.6° for 1019-4 eV, 7.4° 
for 10!9-> eV, and 6.8° for 1019-6 eV. 

The steps for performing the MC simulations are as 
follows. At each energy threshold, 10'94 eV, 10195 eV, 
and 10!°-6 eV, we throw 5x 10° sets of MC trials with the 
same statistics as the data and perform a Li-Ma analysis 
of each trial, which gives us the maximum Li-Ma signif- 
icance and its position. Then, we calculate the angle, 
Omc, between the position that has the MC’s maximum 
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Figure 1. Sky maps in equatorial coordinates using Hammer projections. The color scheme indicates the Li-Ma significance and 
shows the excess (red) or deficit (blue) of events compared to isotropy at each grid point. The positions of maximum excesses 
are marked with the black diamonds. An intermediate angular scale of 20°-radius circles was used for oversampling analysis for 
different energy thresholds. The energy cut for each map is (a) E > 10'°* eV, (b) E > 10*°5 eV, (c) E > 10'*® eV, and (d) 
E > 57 EeV. For three lower energies, (a) to (c), a new excess of events is consistently observed in the same direction. When 
the energy cut is raised to E > 57 EeV (d), the maximum significance moves to the previously observed hotspot, however, a 
smaller excess remains. Note that the right ascension of 0° is at the center of the sky map. 


Li-Ma significance and the center of the PPSC. We count 
as successes those within angle 6p, of the PPSC with 
an equal or greater significance than the data with the 
PPSC: (Sme > Sobs) and (Amc < Oops). By requiring the 
MC set to meet or exceed the conditions of the data, we 
estimate the probability of having an equal or greater ex- 
cess coincident with the PPSC by chance. The chance 
probabilities of having an excess this significant overlap- 
ping the PPSC are estimated to be 3.60, 3.60, and 3.40 
for E > 10194 eV, E > 10°?" eV, and E > 10?" eV, 
respectively. 

We investigate the data further by taking into account 
the major objects in the local large-scale structure of 
the universe similar to the PPSC within the TA’s field 
of view: the Virgo cluster (17 Mpc), Coma superclus- 
ter (90 Mpc), Leo supercluster (185 Mpc), and Hercules 
supercluster (135 Mpc). Their center positions are de- 
fined in the same manner as defining the PPSC center, 
which are marked with the black squares in Figure 4. 


As an example, the Li-Ma significance map for the data 
set with Æ > 109° eV is shown. The data do not show 
an excess at any of the locations of other major objects. 
All Li-Ma significances are less than lo. 

However, to see if it is likely that an excess at one of 
their locations could occur by chance, we repeated the 
MC calculation described in the previous paragraph, re- 
placing the PPSC with all the major structures listed 
here including the PPSC. The chance probability of hav- 
ing an excess of equal or greater significance than that 
of the data on top of any of the five major structures are 
estimated to be 3.10, 3.00, and 2.90 for E > 10194 eV, 
E > 10!°5 eV, and E > 10!%° eV, respectively. These 
significances are sufficiently similar to that of the PPSC 
alone (about 3.50) that we conclude that random coin- 
cidences with the major objects of the local large-scale 
structure are unlikely. 

We repeated the calculation using NASA/IPAC Ex- 
tragalactic Database (NED 2021) centers of the major 
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Figure 2. Sky maps in equatorial coordinates. The Li-Ma significance analysis maps for the data set with E > 10'9° eV. The 
results for the first 5-years of data are shown in the left and for the last 6-years of data on the right. The same color scheme as 
Figure 1 is used. Both maps have a similar, ~3ø, excess in new region which indicates that it may be a steady state. 


Figure 3. Expanded sky maps showing the new excesses of events overlaid with the major clusters and groups of galaxies of 
the Perseus-Pisces supercluster (PPSC). The color scheme is the same as that in Figure 1. The map is shown for three energy 
thresholds with minimum energy increasing from left to right: (a) E > 10*°™4 eV, (b) E > 10° eV, and (c) E > 10'*° eV. 
The representative elements of the PPSC from reference (PPSC 2021) are indicated on the maps with black asterisks. Galaxies 
from the 2MASS Redshift Survey catalog (Huchra et al. 2012), with distances between 35 Mpc and 100 Mpc, are indicated 
with gray dots. These distances are similar to those of the PPSC representative elements. The positions of maximum excesses 
are marked with the cyan diamonds and the center of the PPSC is marked with the blue squares. It is seen that the excess is 
coincident with the overall distribution of the PPSC. The angular separations between the positions of the maximum excesses 
and the center of the PPSC are less than ~10°. 


structures, and the significances were the same. The re- 
sults are summarized in Table 1. Random coincidences 
between the data and the PPSC occur at the ~3.50 
level. This result indicates that it is likely there has a 
cosmic ray source in the Perseus-Pisces supercluster. 


4. SUMMARY 


The TA Collaboration has observed a new excess of 
events in the arrival direction distribution. We found 
the excess over the isotropic background to have local 
significances of 4.40, 4.20, and 4.00 for events of energy 
E > 10"* eV, E > 10 eV, and E > 10°" eV, re- 
spectively, by using the Li-Ma method and a 20°-radius 
circle oversampling analysis. This excess overlaps with 


the Perseus-Pisces supercluster which is a nearby ele- 
ment of the local large-scale structure of the universe 
and is the closest supercluster to us (other than the 
Virgo supercluster within which we reside). 

When looking at the data overlaid with the PPSC, 
the excess is coincident with the overall distribution of 
the clusters and groups of galaxies within the Perseus- 
Pisces supercluster. To determine the probability that 
the data’s 4.00—4.4o significances could occur by chance 
close by the Perseus-Pisces supercluster, we generated 
isotropic Monte Carlo event sets with the same statis- 
tics as the data, thrown according to the acceptance of 
the TA surface detector. The chance probability of the 
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Figure 4. Sky map in equatorial coordinates. The Li-Ma significance analysis map for the data set with E > 10'°° eV. The 
color scheme is the same as that in Figure 1. The representative clusters and groups of galaxies of the Perseus-Pisces supercluster 
(PPSC) are indicated by the black asterisks. The nearby major structures within the Telescope Array field of view—Virgo cluster 
(17 Mpc), PPSC (70 Mpc), Coma supercluster (90 Mpc), Leo supercluster (135 Mpc), and Hercules supercluster (135 Mpc)—are 
indicated by black squares. The data do not show an excess at any of the locations of other nearby major structures other than 


the PPSC. None of them have Li-Ma significances greater than lo. 


Table 1. Summary of the Monte-Carlo studies that estimate the chance probability of having an excess 


Energy (eV) | Events | Criteria | Perseus-Pisces | Any of the Five 
supercluster Major structures 
E > ue 864 4.40 &8.6° 3.60 3.10 
E > 10°53 558 | 4.20&7.4° 3.60 3.00 
E > 10" 335 | 4.00&6.8° 3.40 2.90 


excess of events occurring coincident with the Perseus- 
Pisces supercluster has 3.50 significance. 

We investigated whether there is another excess close 
to the locations of any of the nearby major structures 
similar to the Perseus-Pisces supercluster. None of them 
have Li-Ma significances larger than lo. We repeated 
this process, testing the Monte Carlo event sets against 
the five nearby major astronomical structures. The sig- 
nificance of a random coincidence with any of them is 
estimated to be ~3ø. 

The excess of events observed in the direction of the 
Perseus-Pisces supercluster indicates that a cosmic ray 
source likely exists in that supercluster. The super- 
cluster contains many interesting astronomical objects, 
including active galaxies, starburst galaxies, and large- 
scale shocks, that may be UHECR sources. It is impor- 
tant to study these astronomical objects in the super- 
cluster further, and to increase the statistical power of 
Northern Hemisphere cosmic ray studies. 


The recent TA x4 project, which is the extension of the 
TA SD aperture by a factor of 4 and includes two new 
fluorescence detector stations which overlook the TA x4 
SD array, is designed to study the highest energy cosmic 
rays (Abbasi et al. 2021b). As of 2021, more than half 
of the TAx4 SDs have been deployed and are operating 
successfully. Completing the construction of TAx4 and 
continuing to run TA will likely be an essential key to 
solving the problem of the origin of ultrahigh energy 
cosmic rays. 

A table of events with energies above 1 eV is avail- 
able in a machine-readable form in the online journal. 
It includes date and time, zenith angle, right ascension, 
and declination. 
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